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We investigate the thermal behavior of the (111) surface of silver, using phonon frequencies obtained from ab initio 
total energy calculations, and anharmonic effects treated within a quasiharmonic approximation. Our results 
reproduce the experimental observation of a large and anomalous increase in the surface thermal expansion at 
high temperatures |lj . Surprisingly, we find that this increase can be attributed to a rapid softening of the in- 
plane phonon frequencies, rather than due to the anharmonicity of the out-of-plane surface phonon modes. This 
provides evidence for a new mechanism for the enhancement of surface anharmonicity. A comparison with Al(lll) 
shows that the two surfaces behave quite differently, with no evidence for such anomalous behavior on Al(lll). 



PACS numbers: 68.35, 63.20.Ry, 82.65.Dp 

The equilibrium lattice constant of a crystal is deter- 
mined by a balance between the various attractive and re- 
pulsive forces present within the solid. When the crystal 
is cleaved to form surfaces, this balance is destroyed, and 
the atoms at the surface therefore relax either inwards 
or outwards. Both calculations and experiments show 
that the low-index surfaces of most metals, including the 
noble metals, relax inwards (at temperature T — > K), 
in accordance with simple bond strength - coordination 
number arguments. 

However, recent experiments jj| show that this situa- 
tion is dramatically reversed for the Ag(lll) surface as 
the crystal is heated: Up to T w 670 K, the first inter- 
layer spacing d\% is indeed contracted by ~ 2.5% relative 
to the bulk separation ds', but upon increasing T further, 
d\i increases much more rapidly than ds does, so that 
by 1150 K it is expandedby ~ 10%. Correspondingly, the 
surface coefficient of thermal expansion as becomes more 
than ten times as large as the bulk coefficient of thermal 
expansion as 0- Since thermal expansion arises from 
the anharmonicity of the interatomic potentials, this in 
turn indicates a significant enhancement in anharmonic- 
ity at the surface. 

It has long been realized that it is possible for measures 
of anharmonicity (such as coefficients of thermal expan- 
sion, mean squared displacements of atoms, and the rate 
of change of phonon frequency with temperature) to be 
larger at surfaces than in bulk crystals ||; and the con- 
ventional lore recognizes two possible reasons for such 
an enhancement: (i) the breaking of symmetry due to 
the presence of the surface makes the interlayer poten- 
tial more asymmetric at the surface than in the bulk, 
increasing the size of the odd terms in a Taylor series ex- 
pansion of the energy in powers of atomic displacements; 
(ii) atomic displacements may be greater at the surface 
than in the bulk, thus increasing the relative magnitudes 



of the higher-order anharmonic terms in this series ex- 
pansion. However, early experiments and calculations 
PI , Q on metal surfaces showed that measures of anhar- 
monicity are typically one to three times larger at the 
surface than in the bulk, and thus the significantly larger 
enhancement seen on Ag(lll) is remarkable, and is espe- 
cially unexpected for a close-packed face-centered-cubic 
(fee) (111) surface. 

Moreover, there is a serious disagreement between ex- 
periments and calculations on the Ag(lll) surface: Lewis 
m has carried out molecular-dynamics simulations us- 
ing interatomic potentials obtained using the embedded- 
atom-method (EAM) to investigate the thermal behavior 
of Ag(lll). The results of these simulations differ signif- 
icantly from those reported experimentally: The surface 
layer relaxes inwards at all temperatures, and as is less 
than twice as large as as ■ 

We note that there have been other recent experimen- 
tal observations of large enhancements in surface anhar- 
monicity, but on more open surfaces, viz. Pb(110) [Q, 
Ni(001) H and Cu(110) §. In contrast to the situa- 
tion for Ag(lll), the results of simulations using EAM 
potentials to study Cu(110) and Ni(001) |yj show 
the same qualitative features as the experimental data, 
though the calculated enhancement is somewhat smaller 
than that measured experimentally. 

Why then is there a disagreement between experiment 
and calculation for Ag(lll)? Is it due to inadequacies 
of the EAM potentials? Or could it be an indication of 
some hitherto undetected surface phase transition? 

To study these questions, we have investigated the 
anharmonic properties of Ag(lll) and (for purposes of 
comparison) Al(lll) by performing ab initio calculations 
using density functional theory. Fully separable norm- 
conserving pseudopotentials fl2|| were used in our calcu- 
lations, together with a plane wave basis set with an en- 
ergy cut-off of 60 Ry (20 Ry for Al), and the local-density 
approximation with Ceperley- Alder exchange-correlation 
|l3|. We have verified that these pseudopotentials pro- 
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vide a satisfactory description of harmonic as well as an- 
harmonic properties of bulk Ag and Al . 

The surface calculations were performed using a re- 
peated slab geometry consisting of six atomic layers sep- 
arated by a vacuum layer of the same thickness. The 
k-point sets used to sample reciprocal space contained 
seven points in the irreducible part of the Brillouin zone 
for the undistorted surface; the number of k-points was 
correspondingly increased upon breaking symmetries by 
distorting the lattice in order to calculate phonon fre- 
quencies. Convergence of calculated anharmonic quan- 
tities with respect to energy cut-off, number of k-points 
and number of layers was carefully tested for. 

Our strategy is to compute static energies and phonon 
frequencies by performing self-consistent calculations at 
T = K, and then extend our results to finite tempera- 
tures by using a quasiharmonic approximation. 

In order to obtain a qualitative understanding of the 
mechanisms in operation, and a first estimate of the size 
of surface anharmonic effects, we use a simple model of 
the lattice dynamics of the surface, considering only three 
phonon modes, in all of which the topmost layer of the 
slab moves as a whole - i.e., we assume that the dis- 
placements are confined to the first layer of atoms at the 
surface, and consider only those modes with zero wave- 
vector. We consider one mode in which the surface atoms 
vibrate normal to the surface plane (along the z direc- 
tion), and two modes (along x = [110] and y = [112]) in 
which they vibrate in the plane of the surface. 

We first computed the change in the total energy of the 
Ag(lll) slab upon varying the first two interlayer separa- 
tions di2 and ^3, and maintaining the fee stacking of the 
bulk crystal. This not only provides the static interlayer 
potential, but is also equivalent to simulating the vibra- 
tional mode along z. Our result for the dependence on 
di2 of the first interlayer potential is plotted in Fig. 1; it 
is clearly asymmetric about the minimum at d\2 = 2.30 
A. We found that allowing for the relaxation of c?23 does 
not have a significant impact on the results for the close- 
packed (111) surface |lq| . For the results presented in 
Fig. 1 and the rest of this paper, c?23 is therefore fixed at 
the bulk interlayer separation of 2.34 A. 

To see how this anharmonicity of the interlayer po- 
tential is manifested at finite temperatures, we consider 
a one-dimensional quantum oscillator vibrating in the 
interlayer potential of Fig. 1. A numerical solution of 
the Schrodinger equation for this problem furnishes the 
eigenstates and eigenvalues of such an oscillator, and the 
mean displacement (d\2) n in the n-th eigenstate is ob- 
tained by computing the expectation value of the dis- 
placement operator in each state. The average value at a 
finite temperature T is then obtained by weighting these 
results with the corresponding partition function. 

Our results for d\2(T) obtained from this procedure 
(see the open circles in Fig. 3) indicate a modest en- 
hancement in as relative to cub of ~ 1.7, which is in 



agreement with the EAM results but much smaller 
than that measured experimentally. 

We next performed "frozen-phonon" calculations to 
study the behavior of the two in-plane modes in our 
model. At each value of d±2, the atoms in the surface 
layer were displaced along first the x and then the y di- 
rection, and the total energy was computed for a series 
of displacements up to ±0.15 A. The curvature of the 
resulting plots of energy versus displacement gives the 
mode frequency. We find that the frequency of these in- 
plane modes decreases surprisingly rapidly upon increas- 
ing di2', these results are plotted in Fig. 2. The value of 
the Griineisen parameter 7 extracted from our results for 
the in-plane modes (between 5.5 and 7.5) is significantly 
larger than the average 7 of 2.46 for bulk modes [Elf . 
Thus, the frequently made assumption Q] that 7 is ap- 
proximately equal for surface and bulk modes is clearly 
invalid in this case. 

The surface can thus reduce its vibrational free energy 
significantly by expanding outwards, though such an out- 
ward expansion would be accompanied by an increase in 
the static energy. The optimal value of d\2 is determined 
by minimizing the free energy ^] : 

F(d 12 ,T) = E stat (d 12 ) + ]T F^ b (d 12 ,T), (1) 

i=x,y 1 z 

at each temperature T . Here, E sta ±(di2) is the static 
interlayer potential plotted in Fig. 1, and F* ih (di2,T) is 
the vibrational free energy corresponding to vibrations in 
the i-th direction, which, in the quasiharmonic approxi- 
mation, is given by ||: 

Fi ih (di2,T) = kB Tln{2sinh(^^)}; (2) 

where fcs and h are Boltzmann's constant and Planck's 
constant respectively. The frequency of the mode po- 
larized along the i-th direction, when the first interlayer 
separation is fixed at di 2 , is denoted by uji(d\2); its value 
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FIG. 1. Static interlayer potential between the two outer- 
most layers of Ag(lll). The horizontal lines and open circles 
indicate the energy eigenvalues E n and the mean displace- 
ments (di2)n respectively of an Ag atom vibrating in this 
potential. 
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is obtained directly from our frozen-phonon calculations 
for the two in-plane modes (see Fig. 2). For the out-of- 
plane mode, we compute F* ih (d\2,T) numerically, in such 
a way as to reproduce the exact result for d\2{T) that we 
have already obtained by solving the Schrodinger equa- 
tion when only the mode along z is present. 

Our final result for d\2(T) in the presence of all three 
modes, obtained by minimizing the free energy given by 
Eq. (1), is given by the filled circles in Fig. 3 (note 
that the temperatures are normalized with respect to 
the melting temperature T rn ). The experimentally mea- 
sured data points El are also plotted; both the experi- 
mental and theoretical curves display the same features: 
there is little or no change up to about T/T m = 0.5, i.e., 
T = 617 K, beyond which there is a rapidly increasing 
trend towards outwards relaxation of the surface layer. 
At T/T m = 0.85, i.e., T = 1049 K, we find that the 
surface layer is relaxed outwards by about 15%, whereas 
the experiments show an outwards relaxation of ~ 7.5%. 
Given the simplicity of our model, and the large experi- 
mental error bars, this is as good an agreement as we can 
hope for. 

The increasing slope of d\2{T) reflects a flattening in 
the minimum of the free-energy curve, and the rapid in- 
crease in d\2{T) at high T is a precursor to the devel- 
opment of a saddle-point instability in the free-energy 
curve, similar to that which has been obtained in stud- 
ies of the surface melting of copper surfaces |l7j ]. We 
emphasize that the rapid decrease of u>i{d\2) for the in- 
plane modes is crucial for obtaining the large outwards 
expansion; if, for example, this rate of decrease were to 
be halved, the maximum outwards expansion would be 
drastically reduced to about 2%. 

We expect that our numerical results will change 
slightly upon including other surface phonon modes and 
allowing for dispersion through the surface Brillouin 
zone; however, we note that measurements of surface 
phonon frequency-shifts suggest that the degree of an- 
harmonicity remains approximately constant through a 
surface phonon band llq]. We should also allow for ex- 
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!. Energy fkJi of the top-layer in- plane modes of 
as a function of the first interlayer separation di2- 
Circles and squares indicate modes polarized along the x- and 
y- directions respectively. 
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FIG. 3. Contraction/expansion of di2 relative to ds, as a 
function of normalized temperature T/T m - our calculations 
for Ag(lll) with out-of-plane mode only (open circles), and 
all three modes (filled circles); experiment 1 on Ag(lll) (solid 
line with error bars); our calculation for Al(lll) with all three 
modes (stars). T m is the bulk melting temperature (T^ 9 = 
1234 K, = 933 K). 



pansion in the plane of the surface, but this should in 
fact reinforce the softening of the in-plane modes. 

To check whether such behavior is universal or a pecu- 
liar property of Ag(lll), we repeated the same calcula- 
tions on bulk aluminum and Al(lll). Our results for the 
relaxation of d±2 for Al(lll) are also plotted in Fig. 3, 
and it is obvious that there is no evidence for a dramat- 
ically increased surface anharmonicity on Al(lll). 

Our results indicate that in addition to the two well 
known sources of enhanced surface anharmonicity that 
we have already mentioned, a third (and hitherto ne- 
glected) effect is more important in causing the dramatic 
enhancement in surface anharmonicity on Ag(lll): Not 
only do interlayer potentials at the surface tail off rapidly 
with increasing z , they simultaneously become much flat- 
ter in the xy plane - in other words, the operative effect 
is not so much a decrease in the absolute magnitudes of 
interlayer potentials at the surface, but a reduction in 
their corrugation parallel to surface. As a consequence, 
those surface phonon modes in which atomic dispace- 
ments have significant components in the surface plane 
soften rapidly upon increasing interlayer separations; this 
drives the outermost layer of atoms to expand outwards 
at high temperatures. 

An accurate description of the in-plane corrugation 
clearly requires taking into account the correct distri- 
bution of the electronic charge density (and the result- 
ing chemical bond formation) at the surface. It seems 
plausible that when rebonding effects become significant, 
the EAM (which essentially ignores the relaxation of 
atomic charge densities and the rehybridization of elec- 
tronic states) may fail; this may be why Lewis [[| did not 
observe a large outwards expansion in his simulations. 

The rapid decrease in the corrugation of the interlayer 
potential is evident in Fig. 4, where we have plotted the 
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FIG. 4. Decrease in the corrugation of the interlayer po- 
tential with increasing interlayer separation d\2'- the graphs 
show the increase in surface energy when the outermost layer 
of atoms occupies an atop site or bridge site instead of the fa- 
vored fee hollow site, for Ag(lll) (filled circles) and Al(lll) 
(stars). Arrows indicate the equilibrium value (neglecting 
zero-point vibrations) of d\2 at T — K. 

differences in energies when the outermost layer of atoms 
occupies various stacking sites. Note that (i) the flatten- 
ing occurs more rapidly for Ag(lll) than Al(lll) (ii) 
upon allowing for the lighter mass of Al atoms, the effec- 
tive corrugation relevant for phonon frequencies is actu- 
ally larger for Al(lll) than for Ag(lll). Both these fac- 
tors contribute to the enhancement in 7 of the in-plane 
top-layer modes and thus the larger thermal expansion 
of Ag(lll). 

Such a decrease in the corrugation of the substrate 
potential would also tend to favor a top-layer reconstruc- 
tion of the type that has been observed on Au(lll) Jh| 
or Pt(lll) pc[] ) where the substrate potential is too weak 
to prevent a densification of atoms in the topmost layer. 
However, experiments apparently show no evidence of 
such a reconstruction on Ag(lll) jij; further calculations 
of surface stresses and the strength of intralayer couplings 
should help clarify the situation. 

Our conclusion that the enhancement in surface anhar- 
monicity arises mainly from in-plane vibrations is sup- 
ported by the experimental observation that the magni- 
tude of in-plane vibrations on Ag(lll) rises faster than 
the magnitude of out-of-planc vibrations Exper- 
iments on other surfaces, e.g., Cu(OOl), have detected 
in-plane vibrational amplitudes that are larger than out- 
of-plane amplitudes [ ^2| , and it is interesting to spec- 
ulate whether this (counter-intuitive) result arises from 
the same cause, i.e., from a rapid softening of in-plane 
frequencies due to thermal expansion. 

In conclusion, we have demonstrated that the anoma- 
lously large surface thermal expansion of Ag(lll) can 



be attributed to a rapid softening of in-plane vibrational 
modes (related to a rapid flattening of the corrugation of 
the interlayer potential) upon increasing interlayer dis- 
tances. We have shown that a similar scenario does not, 
however, lead to significant enhancement on Al(lll). 
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